clear;
clc;
a=sym('a%d',[1,4]);
H=hankel(a);
inv(H)